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We study the superconductivity in the three-dimensional multiband d-p model, in which 
a Cu2C>3-ladder layer and a CuC>2-chain layer are alternately stacked, as a model of the 
superconducting spin-chain ladder cuprate. p z -Wave-like triplet superconductivity is found to 
be the most stable, and its dependence on interlayer coupling can explain the superconducting 
transition temperature dependence on pressure in real superconducting spin-chain ladder 
cuprates. The superconductivity may be enhanced if hole transfer from the chain layer to 
the ladder layer can be promoted beyond the typical transfer rate. 
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1. Introduction 

The spin-chain ladder cuprate Sri4_ a; Ca a ;Cu2404i (Srl4-24-41) has attracted much inter- 
est as a related material of high-T c cuprates for more than two decades. It has three types 
of layers, known as the chain, spacer, and ladder. 1-3 ) These layers are stacked in an alternat- 
ing manner along the layer axis b. The chain layer contains Cu02 chains, while the ladder 
layer consists of two-leg CU2O3 ladders. Theoretically, the two-leg spin-1/2 antiferromagnetic 
Heisenberg ladder is a system with a finite energy gap for reaching the lowest magnetic excita- 
tion. 4 ) The hole-doped two-leg spin ladder has been investigated theoretically and confirmed 
numerically to possess an instability towards ci-wave-like superconductivity on the basis of the 
t-J ladder model, 5 ~ 10 ) one-band Hubbard ladder model, 10 ~ 19 ) and three-band Hubbard lad- 
der model. 10 ' 20 ) In 1996, Sri4_ x Ca x Cu2404i with x = 13.6 was found to be superconducting 
under high pressures, P > 3GPa. 21 ) In Sr 14-24-41, doped holes move from the chain layer to 
the ladder layer by Ca substitution for Sr. 22 ) High pressures intensify this self-doping effect 
and afford the electrons in the ladder layer sufficient itinerancy. 23 ' 24 ) However, the supercon- 
ducting transition temperature (T c ) has a maximum at P ~ 5GPa and then decreases with 
increasing pressure. 25 ) The dependence of T c on pressure is common to other compositions as 
well, such as Sri4_. r Ca x Cu2404i for x = 11. 5 26 ) and x = 10, 11, 12. 27 ' 28 ) Such superconductiv- 
ity behavior is similar to that of high-T c cuprates in the overdoped regime. Considering the 
experimental results of the electronic properties, the superconductivity of Srl4-24-41 can be 
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thought of as an extension of that of high-T c cuprates. 29 ) 

In accord with the above novel results, the two-leg Hubbard ladders coupled via a weak 
interladder hopping were investigated theoretically by Kishine and Yonemitsu. 30 ) According 
to their perturbative renormalization group analysis, the system has a <i-wave-like supercon- 
ducting ground state, and restores the interladder coherence within an increase in the extent 
of interladder hopping. Moreover, the superconducting state and its T c were evaluated on 
the basis of the two-dimensional (2D) Trellis-lattice Hubbard model, i.e., the coupled two- leg 
ladder model. 31 ~ 33 ) In their studies, a wide doping regime exists in which d-wave like singlet 
superconductivity appears. In one of them, 32 ) it is also shown that a certain doping regime 
prefers p 2 -wave like triplet superconductivity. Thus, it is plausible that the superconductiv- 
ity may be intrinsically linked to the 2D two-leg ladder structure. Practically, however, the 
superconductivity in Srl4-24-41 appears for certain compounds only under high pressures. In 
order to understand this situation, we should find routes for enhancing superconductivity in 
real compounds, that have not been considered yet in past theoretical studies. It is worth 
noting that Isobe et al. have already pointed out the effects of the hybrid orbital between 
Cu 3d in the ladder layer and O 2p in the chain layer on the superconductivity under high 
pressures. 25, 34 ) According to their arguments, the hybrid orbital accelerates the redistribution 
of holes and enhances the superconductivity. 

In this study, we investigate the superconductivity of the two-leg ladder layer coupled 
with the chain layer. We adopt the three-dimensional (3D) d-p model with the quasi-one- 
dimensional (Q1D) structure in which a CU2O3 ladder layer and a Cu02 chain layer are 
alternately stacked. In our model, the ladder and chain layers are coupled via hybridization 
between the Cu 3d orbital in the ladder layer and the O 2p orbital in the chain layer. Moreover, 
we introduce a sufficiently small on-site Coulomb interaction showing that the second-order 
perturbation theory (SOPT) can be justified. We can treat the superconductivity using a 
weak coupling analysis because, in our model, the effective interaction for Cooper pairing is 
so small that only electrons on the Fermi surface are involved in the superconductivity. In 
particular, the weak coupling formulation by Kondo is applicable even in the case with a very 
small effective interaction. 35 ) We previously applied Kondo's formulation to the study of the 
3D d-p model with multilayer perovskite structure and investigated how the superconducting 
gap depends on the number of layers. 36 ) As in the study of multilayer cuprates, we show that 
calculation on the basis of a 3D model is practical for assessing the superconductivity in spin- 
chain ladder cuprate. Our results give a possible explanation as to why the superconductivity 
appears in Sr^-zCa^C^/tO/u for x > 10 only under high pressures. 

2. Formulation 

Our 3D d-p model is schematically shown in Fig. 1. We consider two Cu sites and three O 
sites in the ladder layer, and one Cu site and two O sites in the chain layer. In the chain layer, 
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the lattice constants along the c-axis and a-axis are \/2 and I/a/2 times those of Sr 14-24-41, 
respectively. Owing to this technique, we do not need a large unit cell, which is needed in 
crystal analysis of Sr 14-24-41, in order to treat the interlayer hopping. We can decompose our 
3D d-p model composed of the ladder layer and chain layer into several parts as follows: 

+ ^Chain + ^Ladder -Chain + H ' 
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where <i/fc CT (^L-) and p^ fccr (p^ fcff ) are the annihilation (creation) operators for d-electrons in 
the Z-th site and for p^-electrons in the m-th site, having a momentum fc and spin a = {t, \}, 
respectively. The site indices I and m, and the orbital index v are defined as shown in Fig. 1. 
In the following, we take both a and b as the unit of length and put a = b = 1. fx represents the 
chemical potential. The noninteracting parts in eq. (1), i.e., H^ ddcr , #g hain , and Ladder-Chain > 
are represented by 
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respectively. In eqs. (2)-(4), we use the abbreviations C| = 2rid p sin4f, = td p e~ lkx / 2 , = 
2ii pp e- ife -/ 2 sin%, = t dp e ifc -/ 2 e fe / 4 , ^ = t dp e ik */ 2 e- ik */ 4 , t% = t pp e ik *, ^ = -t PP e ifc -/ 2 - 
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2t' pp e x cos -f, and 7] k = 2it± sin y-jr ± -j-J- The transfer energies, t^p, t pp , t' pp , and t± are 
denned as shown in Fig. 1. The increase in t± in our model is considered to represent the 
increase in pressure in real Sr 14-24-41. e dp is the level difference between d- and p-electrons. 
Moreover, we use AV to control the charge imbalance between the ladder and chain layers. 
The change in Ay can represent the change in Madelung energy due to the Ca doping of Srl4- 
24-41. Considering only the on-site Coulomb repulsion among d-electrons, the interacting part 
H' in eq. (1) is described as 

u 3 

^' = ^EE d l +q Ak'- q i d ik'i d m- (5) 

1=1 kk'q 

In eq. (5), N is the number of fc-space lattice points in the first Brillouin zone (FBZ). 

In the following analysis, we assume that only electrons on the Fermi surface of the same 
band can have pair instability. For our 3D d-p model, 2 or 3 <i-like bands intersect with the 
Fermi level. Thus, according to the Bardeen-Cooper-Schrieffer (BCS) theory, we have the 
following self-consistent equation for the pair function on the A-th d-like band, 

where i,j = 1,2,3 (Cu site indices) and A, v = 1,2(,3) (d-like band indices). This equation 
is valid for both the spin-singlet and spin-triplet pair functions if we define Vij(q) differently 
as the need arises. Thus, hereafter, we omit the spin indices. Vij(q) represents the effective 
pair scattering between a d-electron on the i-th site and one on the j-th site. The term e v k 
represents the energy dispersion of the i/-th d-like band, and zf(k) represents the matrix 
element of unitary transformation. These variables are obtained by solving the eigen-equation 
for the noninteracting part ^ addcr + tfg hain + ^Ldder- Chain in ^ W- We set = A sc^i 
where A sc denotes the magnitude of and ^ represents its k dependence on the A-th 
d-like band. On the basis of Kondo's argument, 35 ) retaining only the divergent term, we can 
rewrite eq. (6) as 

H = log e A sc • 1 Va(k + k')zt(k)z»{k')5^ kl - (7) 

ijuk' 

for a very small A sc . Equation (7) is a homogeneous integral equation for ^ k with an eigen- 
value of 1/log A sc . We are interested in obtaining the most stable pairing state, so we must 
find the eigenvector ^ k with the smallest eigenvalue l/logA sc using eq. (7) when A sc is 
maximum. Given the quasi-two-dimensionality of our 3D d-p model, i.e., t± <C td p , we assume 
five functions as candidates for the most stable pairing state: 

a x : s-wave-like singlet, (8) 
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In order to solve eq. (7), we substitute these candidates for ^ and and integrate for k z , 
kg . kx t kr£ 5 ky ; Qjiici ky . Then ; we can safely reduce our original eigenvalue problem for ^ to 
an eigenvalue problem for a A [a z (m z ), a x (m z m x ), and a,y(rn y )] in order to obtain the most 
stable pairing state. When we solve it numerically by the standard method, we obtain both 
the eigenvalue l/log e A sc and the eigenvector a x [a z (m z ), a x (m z m x ), and a A (m y )]. Here, to 
solve eq. (7) within SOPT for singlet-pairing states, we have 
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Meanwhile, for triplet-pairing states, 



where 



2T 



(14) 



/£ = 1 [l - tanh 
and T denotes the temperature. 
3. Results and Discussion 

In our present analyses, all e v k and zf(k) in eq. (7) are first calculated for the fc— points on 
an equally spaced mesh in FBZ for each band. The mesh size along k z is 108, and the sizes along 
k^ oLnd ky are both 64. Then, we calculate Vij(k + fc') in eq. (7) only for k— and fc'— points 
satisfying the conditions e k = fi and e u k , = fi, respectively. When we calculate Vij(k + fc') 
according to eqs. (12) and (14), we set the temperature T = 0.001 eV« 10 K, at which our 
system can be considered to behave similarly to the ground state. These calculations have 
been performed at U = 0.3 eV, where magnetic instabilities cannot occur. Other common 
parameters are t dp = 1.00 eV, t pp = -0.50 eV, t' pp = -0.10 eV, and e dp = 2.60 eV. These 
parameters are determined using examples from studies of the three-band Hubbard ladder 
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model, 10 ' 20 ) local-density approximation, 37 ) and angle-resolved photo-emission spectroscopy 
(ARPES). 38 - 40 ) 

As mentioned in the last section, varying the pressure and Ca doping in real Srl4-24-41 
can be reproduced by changing t± and AV, respectively, in our model. Thus, in order to 
comprehensively understand the superconductivity of spin-chain ladder cuprate, we calculate 
the log e A sc of the most stable superconducting state for various t ± and AV values. All through 
these calculations, we fix average holes for Cu 2+ [n h (Cu 2+ )] to 0.25 as well as the value of real 
Srl4-24-41. 1 ) Here, n h (Cu 2+ ) is defined as n h (Cu 2+ ) = [2n£ addcr (Cu 2+ ) + ^ hain (Cu 2+ )]/3, 
where re L a dder(Cu 2+ ) an d n chain(C u2+ ) are average holes per Cu 2+ in the ladder and chain 
layers, respectively. These values should be estimated as the total holes in our model in order 
to be compared with those of real Sr 14- 24-41. ^ We find that the most stable superconducting 
state of our five candidates is always the p 2 -wave-like triplet, represented by eq. (10). We 
confirm that log e A sc estimated for this state with M = 16 does not differ from that with 
M = 15 by more than 8%. The following discussion is therefore restricted to the p z -wave-like 
triplet state with M = 16. 

In Fig. 2, we summarize how log e A sc , n^ addcr (Cu 2+ ) and riQ haiQ (Cu 2+ ) depend on t± for 
various AV values. Figure 2(a) shows the log e A sc dependence on t± for AV = — 0.12eV and 
AV = —0.11 eV, as well as the T c dependence on pressure in real Sr 14-24-41, 25 ~ 28 ) i.e., it has a 
maximum at an intermediate t±. Since the temperature T = 0.001 eV~ 10 K is thought to be 
sufficiently low, T c should be proportional to A SC (T) within the BCS theory. Thus, our results 
of the log e A sc dependence on t± for AV = — 0.12eV and Ay = — O.lleV can qualitatively 
reproduce the T c dependence on pressure in real Srl4-24-41. Moreover, as shown in Figs. 2(b) 
and 2(c), Cu 2+ holes are transferred from the chain layer to the ladder layer with an increase 
in t±, which also qualitatively agrees with the experimental results under high pressures. 25, 34 ) 

On the other hand, for AV = — O.lOeV, log e A sc increases with an increase in t±. For 
AV = — 0.09eV, log e A sc decreases with an increase in t±, but it tends to be much larger 
than those in the other cases. The fact that the behavior of log e A sc differs depending on 
AV can be explained by the configuration of the Fermi surface and the density of state 
(DOS) on it. Let us show the DOS's on the Fermi surface for AV = -0.12, -0.11, -0.10, and 
— 0.09eV in Figs. 3-6, respectively. For AV = —0.12,-0.11, and — O.lOeV, we have two 
quasi-one-dimensional (Q1D) branches of the Fermi surface, which we call "Branch 1" and 
"Branch 2", as shown in Figs. 3-5, respectively. Branch 2 becomes increasingly warped as t± 
increases, and simultaneously the DOS on the branch is enhanced. A large DOS is favorable 
for superconductivity. However, the warped Fermi surface makes Fermi surface nesting worse. 
The effective pair scattering Vij(q), defined by eqs. (12) and (14), is mainly enhanced by 
Fermi surface nesting; thus, the warped Fermi surface suppresses pair instability. Owing to 
these two conflicting effects on the superconductivity, the log e A sc values for AV = —0.12 
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Fig. 2. (Color online) (a) log A sc , (b) < addor (Cu 2+ ), and (c) n h chain (Cu 2 + ) for = -0.12, -0.11, 
-0.10, and -0.09 eV. 



and — O.lleV have maxima at approximately t± ~ 0.035 and 0.030 eV, respectively, when t± 
varies. Although this situation on the Fermi surface is common for Ay = — O.lOeV, the DOS 
on Branch 2 is enhanced more rapidly with t± , as shown in Fig. 5 (note that the color bar scale 
in Fig. 5 is about twice as large as those in Figs. 3 and 4), and such a large DOS surpasses 
the other effect due to the warped Fermi surface in this case. This is the reason why log e A sc 
for AV = — O.lOeV increases with an increase in t±. On the other hand, for AV = — 0.09eV, 
we have three branches of the Fermi surface. Of these, two are Q1D branches and are called 
similarly to those for AV = — 0.12eV. The other branch of the Fermi surface is a quasi-two- 
dimensional (Q2D) branch, newly labeled "Branch 3" in Figs. 6(a) and 6(b). Noting that the 
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Fig. 3. (Color online) Densities of states (DOS's) on Fermi surface for AV = — 0.12eV: (a) for 
tj_ = 0.010 eV, (b) for t ± = 0.030 eV, and (c) for t± = 0.050 eV. 

color bar scale in Fig. 6 is about 100/3 times as large as that in Fig. 5, Branch 3 has an 
extremely large DOS compared with the other branches regardless of t±. This is why log e A sc 
for AV = — 0.09eV always becomes large. However, as shown in Figs. 6(b) and 6(c), Branch 3 
becomes smaller and finally disappears as t± increases. This leads to a decrease in total DOS 
and the degradation of superconductivity. 

When we compare the Fermi surface obtained by our present calculation with that ob- 
served by ARPES, 40) we find that the results for AV = -0.12, -0.11, and -0.10 eV are 
similar. In other words, in real Srl4-24-41, the charge imbalance between the ladder and 
chain is too large to have three branches of the Fermi surface, as observed in our results for 
AV = — 0.09eV. Thus, the superconductivity in spin-chain ladder cuprate may be enhanced 
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Fig. 4. (Color online) DOS's on Fermi surface for AV = -0.11 eV: (a) for t ± = 0.010 eV, (b) for 
t ± = 0.030 eV, and (c) for t x = 0.050 eV. 

if the large charge imbalance between the ladder and the chain can be resolved, for example, 
by varying the elements of the spacer layer without changing their valences. 

Hereafter, we discuss the gap functions of the p 2 -wave-like superconducting state in detail. 
^ on the Fermi surface for AV = — 0.12eV and for AV = — 0.09eV are shown in Figs. 7 
and 8, respectively. For AV = — 0.12eV, ^ changes its sign on the disconnected parts of 
Branch 1, and it does so on those of Branch 2, as shown in Fig. 7. This is similar to the 
triplet state derived by Sasaki et al. 32 ) Furthermore, has no nodes on the connected parts 
of each branch. Thus, this superconducting state is expected to behave as a fully gapped state 
for the microscopically experimental probe. Actually, the relaxation rate Tj - of the NMR 
measurement suggests that the superconducting state of Srl4-24-41 at pressures of 3.5-3.8GPa 
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Fig. 5. (Color online) DOS's on Fermi surface for AV = -0.10 eV: (a) for t x = 0.010 eV, (b) for 
t ± = 0.030 eV, and (c) for t x = 0.050 eV. 

has an s-wave-like character. 41 ' 42 ) Moreover, in their works, the Knight shift of 63 Cu nuclei for 
the ladder derived from high fields shows no change below T c . This fact strongly suggests that 
the superconducting state in real Srl4-24-41 can be a singlet-Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state or a triplet superconducting state, and the latter is consistent with our result 
for Ay = -0.12 eV. 

However, the fully gapped superconductivity for AV = — 0.12eV does not remain for 
AV = — 0.09eV. As shown in Fig. 8, ^ has nodes on Branch 3 or 2. These results indicate 
that our p 2 -wave-like gap function is not robust in terms of the absence of nodes. If the NMR 
measurement can proceed in Srl4-24-41 under different conditions, we may observe T^ 1 for 
the gapless superconductivity. 
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Fig. 6. (Color online) DOS's on Fermi surface for AV = -0.09 eV: (a) for t± = 0.010 eV, (b) for 
t ± = 0.030 eV, and (c) for t x = 0.050 eV. 

Finally, note that the amplitude of ^ on Branch 3 is smaller than that on Branch 1 or 2, 
as shown in Fig. 8. The large DOS on Branch 3 mainly enhances the superconductivity on the 
other branches, not on itself. This synergistic effect is caused by the interband interaction, 
i.e., Vij(k + k')z^(k)Zj(k') for A ^ v in eq. (6), originating from the mixing between d- 
and p-orbitals. Thus, the hybridization effect also plays a significant role in enhancing the 
superconductivity in this material. 

4. Conclusions 

We have demonstrated that the 3D d-p model with the Q1D structure similar to Srl4-24-41 
can represent p 2 -wave-like triplet superconductivity up to the second order in the perturbation 
theory framework. On the basis of this model, we can reproduce the T c dependence on pressure 



12/16 



J. Phys. Soc. Jpn. Full Paper 




Fig. 7. (Color online) on the Fermi surface for Ay = -0.12 eV: (a) for t± = 0.010 eV, (b) for 
tj_ = 0.030 eV, and (c) for t± = 0.050 eV. 



by changing the interlayer coupling. The calculated results on the Fermi surface configura- 
tion and the superconducting state can give a comprehensive picture to explain the ARPES 
and NMFl experimental results. Moreover, our results show the possibility of enhancing the 
superconductivity if the charge imbalance between the ladder and the chain can be decreased 
by varying the elements of the spacer layer. This speculation is based on the fact that the 
hybridization effect due to interlayer coupling and the other mixing between d- and p-orbitals 
is crucial to enhance the superconductivity. 
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